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Рассматривается возможность применения биоинспирированных методов для решения задачи криптоанализа 
асимметричных алгоритмов шифрования на основе факторизации составных чисел. Представлены алгорит- 
мы муравьиных и пчелиных колоний для разложения составных чисел на множители путём определения 
Делителя числа с заданной точностью в заданном интервале. Описаны отличительные особенности пред- 
ставленных методов, в том числе возможность эффективной параллельной реализации. 
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Введение. В настоящее время в науке и технике находят широкое применение алгоритмы, осно- 
ванные на природных системах. Это новое научное направление под названием «природные вы- 
числения» объединяет математические методы, в которых заложен принцип природных механиз- 
мов принятия решений. К ним относятся методы моделирования отжига, генетические, эволюци- 
онные методы, алгоритмы роевого интеллекта и др. [1]. 

Однако, как отмечено в [2], недостатком эволюционных методов является использование 

«слепого» поиска, что в общем случае приводит к следующим проблемам: генерация решений с 
нарушениями и, как следствие, увеличение времени поиска и необходимость дополнительного 
контроля; генерация большого количества одинаковых решений; генерация большого количества 
плохо приспособленных решений, что в общем случае может привести к попаданию в локальный 
оптимум. Поэтому представляет интерес применение эвристических методов, инспирированных 
природными системами, в которых осуществляется поэтапное построение решения задачи (т. е. 
добавление нового оптимального частичного решения к уже построенному частичному оптималь- 
ному решению). К методам данного вида относят муравьиные и пчелиные алгоритмы, основные 
идеи и принципы которых описаны, например, в [1]. Известны случаи применения алгоритмов 
роевого интеллекта для оптимизации широкого круга задач, в том числе в криптоанализе. Вместе 
с тем следует заметить, что наряду с классическими симметричными алгоритмами шифрования в 
настоящее время используются методы асимметричной криптографии — сравнительно молодой 
области науки. Один из её первых алгоритмов — КВ$А [3, 4], сложность которого определяется 
трудностью факторизации больших чисел. Как отмечено в [5], факторизация (задача разложения 
числа на простые множители) — одна из основных теоретико-числовых задач, используемых в 
криптографии. 
Постановка задачи криптоанализа асимметричных криптосистем. Как отмечено в [6], во 
всех асимметричных криптосистемах (в т. ч. К$А) используются два ключа: один для шифрования 
(Кок), другой для дешифрования (К). Ключи представляют собой пары (К, М), (Кек, №), где 
№ — модуль, при этом М = Р: О, где Ри О — случайные большие простые числа. При этом воз- 
можны несколько вариантов криптоанализа: 
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1. На основе известного открытого текста и шифртекста подобрать такой секретный ключ 
Кекр, Чтобы выполнялось равенство М, = С*(тод №), где М‚— открытый текст, С; — шифртекст. 

2. Опытным путём определить функцию Эйлера $ (Л) = (Р- 1)(О0 - 1), разложив модуль № 
на множители Ри О, и секретный ключ Хек ИЗ СООТНОШения Кью ° Кок = 1(то4 ф (Л). 

3. Подобрать такие числа Ри О, чтобы выполнялось соотношение М = Р` О. 

Поскольку в варианте 1 нахождение секретного ключа К имеет комбинаторную слож- 
ность (в общем случае может потребоваться полный перебор), как и в варианте 2 (требуется пол- 
ный перебор всех взаимно простых чисел в интервале [1, ф (Л/], то актуальна задача исследова- 
ния возможности применения современных технологий природных алгоритмов для определения 
всех возможных делителей составного числа М. 

Отметим, что ранее в [7] был представлен генетический алгоритм (ГА) для решения зада- 
чи определения вариантов разложения заданного числа /М на множители. Однако при реализации 
данного алгоритма основной также является задача определения числа Иа/ являющегося делите- 
лем числа М. Для решения этой задачи в [7, 8] предлагается алгоритм, сущность которого заклю- 
чается в определении на отрезке [0, Л методом генетического поиска точки А, удовлетворяющей 
условиям (№ - А) / К = л— целое (т. е. (М- А) / К- [(М- В) / К =0), 

В: (п+1) =№М, (1) 

где [х] — целая часть числа х, или множества точек К, А, А, .., Ак удовлетворяющих условиям 

К /[ В = т, (№ - К) [| К = №» (® - К») [ К = 1, ..., (Кк- К) | ВЕ = т Вх 

х (++... + пу = М где К, пк — целые числа. Очевидно, что выражение (1) может быть 
использовано в качестве целевой функции в ГА. 

Однако основной проблемой при реализации данного ГА является нахождение экстремума 
немонотонной функции, т. е. функции, значение /(х) которой в каждой точке х является, по сути, 
случайной величиной и не даёт информации о приближении к глобальному экстремуму. Таким 
образом, при нахождении экстремума функции Ё(Ю) = (М- К) / В - [(№М- К) / Ж]| (1. е. значения 
Ю, для которого Ё(Ю) = ((М- К) / К - (М- №) / К] = 0)) с использованием структуры ГА фактиче- 
ски имеет место «слепой» поиск, и, как отмечено в [2], это является основным недостатком 
большинства структур ГА. 

На основе логарифмического закона распределения простых чисел [9] в работе [6] описа- 
на модель ГА для нахождения разложения составного числа М№на 2 простых множителя. В данной 
модели на отрезке [0, ЛМ осуществляется генерация популяции простых чисел. Случайным выбо- 
ром битов генерируется случайное п-битовое число С. Далее на отрезке [6 — п (С), б+п (6])] 
осуществляется поиск наиболее вероятного простого числа. Каждое х е [б -— п (6), бС+ п (6)] 
последовательно проверяется на делимость с простыми числами в интервале [3, 2 ` п (6)]. Так 
как реальное значение С составляет порядка С = 2°1?, то радиус поиска представляется как г = 
= п Сб= 1092С / |0д2е = 09»б / 1,442695 = п/ 1,442695. 

Далее к полученной популяции простых чисел применяется множество генетических опе- 
раций для получения простых чисел-потомков. Для этого используется описанная методика поис- 
ка наиболее вероятного простого числа в окрестности декодированного числа. Однако, учитывая 
реальную размерность задачи (как отмечено в [3, 6], разработчикам криптоалгоритмов на базе 
В5А приходится применять числа длиной не менее 200 десятичных знаков, на данный момент 
криптостойким считается ключ размерностью 217“), этот алгоритм может потребовать значитель- 
ных временных ресурсов (при операции формирования большой популяции простых чисел и их 
получении при проведении генетических операций) и оказаться достаточно трудоёмким при прак- 
тической реализации. Отметим также, что при использовании известных технологий распаралле- 
ливания ГА (например, «островной» модели), описанных, в частности, в [8, 10], необходимо орга- 
низовать межпроцессорные связи между «островами» (группами процессоров, моделирующих 


1545 


Физико-математические науки 








развитие популяции), чтобы предотвратить потери хороших решений и попадание в локальный 
оптимум. Как отмечено в [10], частота миграции является наиболее существенным фактором и 
должна быть установлена на основе экспериментальных результатов, что также может потребо- 
вать временных ресурсов и снизить эффективность работы ГА. 

Поэтому представляет несомненную актуальность разработка новых методов, использую- 
щих модели процессов живой природы и ликвидирующих отмеченные недостатки ГА путём отме- 
ченного выше поэтапного построения оптимальных решений. Рассмотрим возможный подход для 
решения задачи разложения составного числа на простые сомножители с помощью другого клас- 
са биоинспирированных методов — алгоритмов муравьиных и пчелиных колоний. 
Разработка метода муравьиных колоний для факторизации составных чисел. Основные 
положения теории муравьиных алгоритмов и описание их работы приводятся в [1, 11, 12, 13]. В 
соответствии с [1], для того, чтобы построить подходящий муравьиный алгоритм для решения 
задачи, необходимо представить задачу как набор компонент и переходов или как набор неори- 
ентированных взвешенных графов, на которых муравьи могут строить решения. Поэтому рассмот- 
рим возможный подход для сведения задачи разложения составного числа на множители к задаче 
нахождения кратчайшего пути в графе. Пусть задано достаточно большое число Л, которое необ- 
ходимо проверить на простоту и определить его делители на отрезке [п, п] с заданной степенью 
точности, т. е. найти такие числа х, для которых ((/М№/ х) — [№/ х]) >> тт. Пусть числа хк е [п, п 
содержатся в оперативной памяти одного процессора (или «острова» процессоров) и являются 
вершинами полного графа С = (Х И), где |М = п - + 1 — множество вершин, | = 
= (п-т + 1) (11 - п) / 2 — множество рёбер, при этом вес вершины х, равен 2 (х,) = (№М/ х) — 
— [№М/ х], т. е. равен дробной части частного, полученного от деления М№на х,. Сформулируем за- 
дачу нахождения кратчайшего пути в следующей форме: найти маршрут Т в графе С, содержа- 
щий заданное число л7 вершин, удовлетворяющий условию 

У Е(х,) > тп 
1=1,2,..т 
Вес каждого маршрута 7х определим как ([„, РЁ), где 
Е: УЕ) Е, = т (Е (х,)). 


Хх, ЕТ, 
1: К 
Х,ЕТь 


Таким образом, задачей муравьиного алгоритма в данном случае является определение 
маршрута, имеющего минимальный вес /х и содержащего 7 вершин и вершину графа хх, для ко- 
торой Е (жк) = (М/ жд — [М/ ж > тип, т. е. которая является наиболее точным делителем числа 
№. Так как цикл по времени жизни колонии продолжается до тех пор, пока не будут просмотрены 
все вершины графа С, то вершина хх е [п, п] будет определена за число временных циклов С 
алгоритма не менее 

П, -П, +1 
ОЕ Ы (2) 
т.М 
где М — размер популяции муравьёв, л7 — число вершин в маршруте. 

Предположим, что муравьи обладают свойствами поведения (память, зрение, обоняние), 
описанными в [11, 12]. Так как веса вершин Е (х,) априорно неизвестны и определяются в процес- 
се работы алгоритма, то вероятность перехода муравья К из вершины х; в х, устанавливается на 
основе соотношения 


Р(-< (3) 


где к — множество вершин, доступных муравью К, находящемуся в вершине / т, (8) — уровень 
феромона в момент времени {на ребре хх, Е— параметр цикла времени жизни колонии. 
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После прохождения А-м муравьём маршрута из 7 вершин подсчитывается длина пути. Она 
равна сумме всех весов вершин, по которым прошёл муравей. На каждом ребре хх; маршрута ко- 
личество отложенного феромона составляет 


Ат (#) =@/Ь (6), 
где О — параметр порядка длины оптимального пути (определяющий уменьшение Ат с увели- 


чением длины маршрута /х (6). 
Обновление феромона производится в соответствии с выражением 


ь (+) -@-р).[ (= Ульи (.) (4) 


где Е — множество муравьёв, использовавших в маршруте ребро хх, р — интенсивность испаре- 
ния (в [13] предлагается выбирать р = 0,6). 

Отметим, что в соответствии с [12] на начальном этапе задаётся начальное расположение 
муравьиной колонии. В общем случае для данной задачи может быть использована стратегия 
«дробовика», когда ограниченное количество агентов М случайным образом размещается в вер- 
шинах графа без повторений. При такой стратегии в оптимальном случае число циклов жизни 
колонии Сопределяется выражением (2). 

Итак, муравьиный алгоритм для задачи нахождения делителей числа ММ можно сформули- 
ровать следующим образом. 

1. Задание начальных значений параметров О, М, р, вк. 

2. Определение минимальных весов рёбер на основе начальной концентрации феромона. 

3. Задание оценки /” — длина кратчайшего маршрута, оценки Р” — точность решения. 

4. Цикл по времени жизни колонии: Ё= 1. 

5. Размещение популяции М муравьёв в случайные вершины графа без совпадений. 

6. Цикл по муравьям: К = 1. 

7. Муравей К строит маршрут Тк (8) длиной [к (8 на основе распределения вероятности по 
рёбрам в соответствии с формулой (3). 

8. К = К- 1, если К < М, переход к 7. 

9. Проверка всех маршрутов 7» на лучшее решение по сравнению с Гир. 

10. В случае, если получен маршрут Т» (8) длиной /х (8), которая предпочтительнее [^, и 
маршрут содержит вершину х„ для которой Е (х) < Е’, обновить [и Р°. 

11. Цикл по рёбрам графа / = 1. 

12. 1= 1. 

13. Обновить неминимальные следы феромона на ребре х, в соответствии с (4). 

14. } = ]+ 1, если ]< п;- п;+ 1, перейти к 13. 

15. /= 7+ 1, если /< п;- п, перейти к 12. 

16. Е= Ё+ 1, если условия остановки не выполнены, переход к 5. 

17. Вывести маршрут кратчайшей длины [х, а также вершину х, для которой Е (х) мини- 
мально. 

Условиями остановки могут быть Ё> &„«, прекращение минимизации длины маршрутов, 
определение вершины х, для которой РЁ (х) = 0. 

Таким образом, данный алгоритм позволяет определять маршрут, содержащий заданное 
число вершин х, для которых >». Е (х) = тт (т. е. являющихся самыми точными делителями числа 
М№Мна отрезке [л, п]), а также найти вершину х, для которой Е (х) = п, т. е. которая является 
делителем числа №с заданной степенью точности. 

Рассмотрим пример. Пусть /М = 893 и в памяти процессора (или «острова» процессоров) 
содержится отрезок [5, 45], т. е. в этом случае |Х| = 41. Определим М = 4 (количество муравьёв), 
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т = 4 (число вершин в маршруте), О = 4. Начальную концентрацию феромона при ЕЁ = 0 зададим 
равной т, (0) = 0,25. Для решения задачи нахождения маршрута, содержащего л7 вершин, и зна- 
чения х, для которого Р (х) = тт, разместим М муравьёв в случайно выбранные вершины графа, 
например в вершины 5, 10, 30, 39. Так как на Ё = 1 итерации веса рёбер одинаковы, определим 
случайным образом 4 маршрута. 


П 4): 5—12.—18—25 1 (1) = 0,6 + 0,42 + 0,61 + 0,72 = 2,35 
ТЬ (1): 10—22 —18 — 36 5 (1) = 0,3 + 0,59 + 0,61 + 0,81 = 2,31 
7: (1): 30 —25— 17—40 13 (1) = 0,77 + 0,72 + 0,53 + 0,33 = 2,35 
(1): 3937—1124 1 (1) = 0,90 + 0,14 + 0,18 + 0,21 = 1,43 


Таким образом, на рёбрах маршрута К количество отложенного феромона после 1 итера- 
ции ти (0) + Ах = ти (0) + 9/4, (1) составит: 

и (0) -- Ат; =0,25+1,7 =1,95; ет (0) -- Ат; =0,25 +1,73 =1,98; 
тя (0) + ху =0,25+1,7=1,95; ти (0) + Ах; =0,25+2,8 = 3,05. 

После испарения получим следующие значения концентрации на рёбрах графа, соответ- 
ствующих маршрутам: т, (1) =0,78; т, (1)=0,79; т (1)=0,78; ту (1)=1,22. Рёбра графа после 
1-й итерации будут иметь веса: т512 (1) = тизлв (1) = т1в5 (1) = тзо5 (1) = 12547 (1) = т17ло (1) = 
= 0,78; 110,22 (1) = 122,18 (1) = 118,36 (1) = 0,79; 739,37 (1) = 137/11 (1) = 11124 (1) = 1,22. При ЭТОМ 
Е = тю (Е (х,)) = 2 (12)=0,42; = Е (10) = 03; В = Е (40) = 0,33; Е, = Е (37) = 0,14; Р' = 


Х, ЕТ, (1) 
= 0,14. Для остальных рёбер веса т, (1) = 0,25. 
Перейдём к следующей итерации при Ё = 2. Разместим муравьёв в случайно выбранные 
позиции, например, в 33, 42, 44, 10, и пусть определены следующие маршруты. 


Т; (2): 33 — 27 — 37 — 11 Ё (2) = 0,06 + 0,07 + 0,14 + 0,18 = 0,45 

Т> (2): 42 —9 — 39 — 37 Ь (2) = 0,26 + 0,22 + 0,90 + 0,14 = 1,52 

Тз (2): 44 — 10—26 —27 [3 (2) = 0,30 + 0,30 + 0,35 + 0,07 = 1,02 

Та (2): 10 — 13—25 —17 [4 (2) = 0,30 + 0,69 + 0,72 + 0,53 = 2,24 
Таким образом, на рёбрах маршрутов количество отложенного феромона составит: 

Алу = 9/1 (2) = 4/0, 45 = 8, 89; Ата = О /Ь (2) =4/1,52 =2,63; 

Ату = О/Ь (2) =4/1,02 = 3,92; Ати = 0/1, (2) = 4/2,24 =1,79. 


Таким образом, после испарения феромона рёбра графа будут иметь следующие веса: 
15,12 (2) — 11218 (2) — 118,25 (2) — 130,25 (2) = 125,17 (2) — 117,40 (2) = 0,31; 110,22 (2) — 12218 (2) > 
= 11836 (2) = 0,32; т1олз (2) = 113,25 (2) = 12517 (2) = 0,82; таало (2) = т1орв (2) = 126,27 (2) = 1,67; 
1429 (2) = тозо (2) = 1,15; тзэз7 (2) = 1,54; тззр7 (2) = т27з7 (2) = 3,66; тзу11 (2) = 4,04; тим (2) = 
= 0,49. 

После 2-й итерации А, = Ё (33) = 0,06; Р = Е (9) = 0,22; В = Е (27) = 0,07; А = Е (10) = 
= 0,3, Е’ = 0,06. Для остальных рёбер веса ту (2) = 0,25. 

Таким образом, данная итерация показывает, как маршрут, проходящий по вершинам с 
наименьшим весом, максимально обогащается феромоном. 

Перейдём к итерации Ё = 3 и разместим муравьёв в вершинах 19, 8, 32, 10. Пусть опреде- 
лены следующие маршруты. 


РЗ -- 97 —3 1; (3) = 0,0 + 0,14 + 0,07 + 0,06 = 0,27 

т (3):8=—13=:40—22 1. (3) = 0,63 + 0,69 + 0,30 + 0,59 = 2,21 
7; (3): 32 —33—27— 37 13 (3) = 0,91 + 0,06 + 0,07 + 0,14 = 1,18 
№3): == 937 1 (3) = 0,30 + 0,53 + 0,90 + 0,14 = 1,87 
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Таким образом, после 3-й итерации количество отложенного феромона на рёбрах, входя- 
щих в маршруты, составит 


Ату = @/ (3) = 4/0,27 =14,81; ту = О/Ь (3) =4/2,21 =1,80; 
Ата = О /Ь (3) = 4/1,18 = 3,39; Ати = 9/1, (3) = 4/1,87 = 2,14. 


Таким образом, после испарения феромона рёбра графа будут иметь следующие веса: 
15.12 (3) = т12лз (3) = 11825 (3) = тз0,25 (3) = 125117 (3) = т174о (3) = 0,124; т1ор2 (3) = 0,84; т›2 лв (3) = 
= 118,36 (3) = 0,128; 739,37 (3) = 1,47; 137/11 (3) = 1,62; 211,24 (3) = 0,20; 742,9 (3) = 19,39 (3) = 0,46; 
14410 (3) = 110, (3) = 12627 (3) = 0,67; т1з,25 (3) = 12547 (3) = 0,33; тлэзу (3) = 6,02; тз727 (3) = 
12733 (3) = 8,74; 78,13 (3) = 0,82; 713,10 (3) = 1,05; 132,33 (3) = 1,46; 710,17 (3) = 51739 (3) = 0,96. 

После 3-й итерации А; = РЁ (19) = 0; Б = Е (10) = 0,30; А = Е (33) = 0,06; Е = Е (37) 
0,14; Е” = 0. Для остальных рёбер веса х, (3) = 0,25. 

Таким образом, данный пример демонстрирует, как на 3-й итерации маршрут, содержа- 
щий вершины с наименьшим весом (19 — 37 — 27 — 33), наиболее сильно обогащается феромо- 
ном. Значения, соответствующие этим вершинам, являются, очевидным образом, делителями 
числа Л с максимальной степенью точности = (для данных вершин 0 < = < 0,14). Очевидно, что на 
последующих итерациях по этому маршруту пойдёт наибольшее количество агентов, в то время 
как другие пути будут исчезать. 

Таким образом, отличительными особенностями данного подхода являются: 

- возможность эффективной параллельной реализации, связанной с отсутствием миграции 
особей между процессорами (в отличие от классического ГА); 

- возможность определения множества значений х, являющихся наилучшим приближени- 
ем делителя числа //с заданной степенью точности. 

Отметим, что, как следует из приведённого примера, веса вершин графа могут опреде- 
ляться в процессе работы алгоритма (при формировании маршрута муравьёв), т. е. после этого 
задача сводится к реализации классического муравьиного алгоритма для определения кратчай- 
шего маршрута заданной длины. 

Экспериментальные результаты показывают, что эффективность муравьиных алгоритмов 

растёт с увеличением размерности решаемых задач оптимизации. При этом сходимость и качест- 
во решения зависят от начального расположения колонии и выбранных параметров. В качестве 
теоретической оценки размерности колонии для данной задачи можно принять 
М> (п; — п + 1) / т. Так как в этом случае маршруты муравьёв с высокой степенью вероятности 
пересекаются, также высока и вероятность того, что муравей повернёт на оптимальный путь, 
пройдёт по нему и обогатит его феромоном. 
Алгоритм разложения составных чисел на простые сомножители с использованием 
пчелиных колоний. Необходимо решить следующую задачу. На отрезке [п, п] определить це- 
лочисленные делители числа // являющиеся простыми числами, т. е. осуществить разложение 
числа /Ина простые множители. Поскольку простые числа, как отмечено в [6, 9], на заданном ин- 
тервале распределены по логарифмическому закону, то целью поиска является определение в /- 
окрестности точки х всех простых чисел у; и точности =, с которой данное простое число у; явля- 
ется делителем числа Л. Как и в предыдущем случае, значение = определим как значение функ- 
ции Р (у) = (М/ У) — [М/ У. 

Рассмотрим применение для реализации данного подхода алгоритмов пчелиных колоний, 
являющихся относительно новым «природным» алгоритмом, используемым, как отмечено в [14], 
в первую очередь, для оптимизации сложных многомерных функций. Так как в данном случае 
ищется экстремум немонотонной функции 2 (х), поэтому исследование возможности применения 
для решения данной задачи эвристических методов, не использующих непосредственным образом 
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аппарат математического анализа, является, несомненно, актуальной задачей. Отметим, что опи- 
сание алгоритма, основанного на поведении колонии пчёл, приводится в [14, 15, 16]. Исследова- 
ние пчелиных алгоритмов для решения комбинаторных теоретико-графовых задач (задача раз- 
биения графа, раскраска графа, сравнение с другими биоинспирированными методами) приво- 
дится в [17, 18]. Таким образом, на основе математической модели алгоритма, основанного на 
поведении колонии пчёл, и его описания в [16] алгоритм факторизации числа сформулируем в 
следующей форме. Как и ранее, будем предполагать, что поиск простого делителя х; осуществля- 
ется на заданном отрезке [п, п). 

1. Определить параметры алгоритма: количество пчёл-разведчиков 0, количество рабочих 
пчёл В, количество участков для исследования окрестностей 2, точность нахождения делителя с. 

2. Выбрать на отрезке [п, п] Р значений аргумента хи, ..., Хо. 

3. В выбранные точки х; направить В рабочих пчёл для поиска в их г-окрестности простых 
чисел в соответствии со следующим алгоритмом [6]. 

3.1. Определить для каждого значения х, значение окрестности г = л / 1,442695 (где п — 
число бит в двоичной записи числа). 

3.2. Каждое число у е [х, — г, х› + П последовательно проверяется на делимость с просты- 
ми числами в интервале [2, 2 * 1. 

3.3. К числам, которые прошли тест проверки делимости, применяются известные тесты 
проверки простоты числа, описанные в [19], например, тест Миллера — Рабина, алгоритм, осно- 
ванный на матрице Сандарама [20]. 

4. После определения множества простых чисел У для каждого у; е У определить значе- 
ние функции Ё (у), найти пир А (у,). Определить у, для которых Ё (у) < =. 


5. Из множества У выбрать случайным образом 2 элементов, данные значения обозначить 
как хи, ..., хх. Отправить О пчёл-разведчиков для поиска на отрезке [п, п] О значений аргумента 
Ха, ..., Ха+р. ЕСЛИ условия остановки не выполнены, переход к 3, иначе к 6. 

6. Конец работы алгоритма. 

Условиями остановки алгоритма могут являться окончание временного ресурса; опреде- 
ление величины (или множества величин) х, для которых РЁ (х) = 0 или РЁ (Х) < =; определение 
значений функции 2 (х) для всех х; е [п, п. 

Таким образом, в данном алгоритме выбор значений аргумента х; = [л, п] имитирует по- 
ведение пчёл-разведчиков, а поиск в /-окрестности наиболее вероятных простых чисел имитирует 
поведение рабочих пчёл (пчёл-фуражиров). Поскольку в данном случае определяется экстремум 
немонотонной функции, то выбор точек х, на отрезке [лп, п] для поиска простых чисел в их окре- 
стности производится на каждой итерации случайным образом, что приводит в общем случае к 
равновероятной возможности получения глобального оптимума на каждой итерации (в отличие 
от направленного схождения к экстремуму в классическом пчелином алгоритме, описанном в 
[16]). 

Отличительной особенностью разработанного алгоритма, как и в [16], является динамиче- 
ское разбиение поискового пространства на случайные области, что, с одной стороны, уменьшает 
время работы алгоритма, а с другой — повышает вероятность нахождения глобального оптимума 
на каждой итерации. В то же время благодаря случайному равновероятному поиску по всей длине 
отрезка становится возможным применение эффективных стратегий распараллеливания, приво- 
дящих к сокращению временных затрат. На каждой итерации временные затраты равны поиску в 
самой большой /-окрестности. 

Рассмотрим пример применения итерационного алгоритма. Пусть № = 15589 и в памяти 
процессора содержится отрезок [п, п] = [115, 140]. Определим Д = 2 и выберем случайным об- 
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разом точки х, = 123, ж = 133. Так как ж = (123) = (1111011)>, то определим п = 7иг= 
= п/ 1,442695 = 5. 

Выделим интервал чисел [118, 128]. Проверим делимость чисел в данном интервале на 
простые числа в интервале [2, 2 ` /| = [2, 10], т. е. на {2, 3, 5, 7}. Получим числа {121, 127}. Од- 
нако, т. к. число (121 - 1) / 2 = 60 содержится в матрице Сандарама в позиции (5, 5), то число 
121 является составным. Для числа 127 Е (127) = 0,75. 

Так как № = (133) = (10000101)›, то п =8иг= 6. 

Рассмотрим интервал чисел [127, 139]. Проверим делимость чисел из интервала на про- 
стые числа в интервале [2, 12], т. е. на {2, 3, 5, 7, 11%. Получим числа У = {127, 131, 137, 139%. 
Так как для всех уе У(у- 1) / 2 не содержится в матрице Сандарама, то данные числа являются 
простыми. При этом Р (127) = 0,75; Е (131) = 0,0; 2 (137) = 0,79; Е (139) = 0,15. Таким образом, 
на данной итерации определено значение х; е [п, п], являющееся точным делителем числа М. 

Отметим, что в общем случае задачу проверки простоты числа ЛМ с помощью матрицы 
Сандарама можно свести к проверке принадлежности числа № = (М- 1) / 2 множеству арифме- 
тических прогрессий, составляющих строки (столбцы) матрицы, приведённой в [20]. То есть зада- 
ча сводится к определению целочисленного л из соотношений 
4+3 (п- 1) =м№7+5'(п-1=М№М..., ак+ @` (п- 1) = №, где а=4+3` (К- 1), Ч=З+ 
+2: (К- 1), т.е. п- 1= (№- (4+3: (К- 1))) / (3+2: (К- 1). 

Поскольку определение целого значения л имеет смысл при 

№ - (4+3. (к-1)) т 
3+2. (К-1) 
то определим значение К — число прогрессий (строк или столбцов матрицы), которые надо про- 
верить для определения, является ли № основой простого числа (в предположении, что № = 4 не 
является первым членом первой прогрессии). Преобразуя (5), получим № - 3 К-1>2'К+1, 


откуда 
№'-2 
К Е 6 
< - (6) 


где ]х{ — ближайшее целое снизу к числу х. Таким образом, если число № не является членом 
первых К арифметических прогрессий матрицы Сандарама (где К определяется выражением 
(6)), то число М достоверно является простым. Отметим, что для предлагаемого в [6, 19] ве- 
роятностного теста Миллера — Рабина рекомендуется в общем случае проверка д = 1092 
свидетелей простоты, после чего вероятность того, что М составное, не превышает 4“. 
Поскольку наиболее целесообразными являются криптосистемы, в которых простые со- 
множители Ри О, составляющие модуль М, имеют порядок 2°1?, то максимальное значение г = 
= 355, поэтому длина интервала |2, 2 `/| = 700. Простые числа в этом промежутке могут быть 
найдены с помощью известных таблиц, поэтому их определение не влияет на временную слож- 
ность алгоритма. 
Выводы. Таким образом, в данной работе был представлен возможный подход к факторизации 
составных чисел с использованием биоинспирированных методов (алгоритмов муравьиных и пче- 
линых колоний), описаны основные отличительные особенности методов (в частности, возмож- 
ность эффективной параллельной реализации), были также представлены демонстрационные 
примеры, иллюстрирующие возможность практического использования предложенных алгорит- 
мов. 


(5) 
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СОМРОЗТТЕ МУМВЕВ РАСТОКТ2АТТОМ 


А. $. бегдеуеуи 

(Роп Зае Теспткса! Упмег$КУ), 

О. Р. Тгебуакох, А. Е. УазПуеуи 

(Кгазподаг Бгапсй оЁ МИЙКагу Асадету оЁ СоттитсаНоп), 
у. О. СпегпузВем 

(Роп Зае Тесптка! Упмег$Ку) 


Тре арр/саНоп оЁ Те Мотзр!иеа те тоа$ юг вапа!тда те сгурапа$5 ргоМет оЁ Те азуттейс епсгурНоп 
а/дотйт5 оп {Пе Ба$/5 о {те сотрозйе питрег гасюпгаНоп 15 соп/аегеа. Тре а/боттб о! апЕ апа Бее соотез 
Юг ше сотрозйе питВЕг Габюпгайоп Бу {пе аейтоп оЕ {пе 1щедег амбог {0 пе зресйеа ассигасу т {пе чаеа 
Иегуа! аге аадисеа. Тве ргорегиез оЁ {те те тоа$ ргезещеа, псиата {те ейсепЕ рага!Е! Геа5 Пу, аге ае- 
$СПБЕД. 

Кеуигога&: сгур{апа/уз5, Бее а/дотё т, апЁ соопу адотт, рвеготопе, асюопгайоп оЁ питрегз, азутте с 
сгурюзу&ет5, Мотзриеа тетосб. 
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